Fractal formation and ordering in random sequential adsorption 
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We reveal the fractal nature of patterns arising in random sequential adsorption of particles with 
continuum power-law size distribution, P{R) ~ R"~^ , R < -Rmax- We find that the patterns become 
more and more ordered as a increases, and that the Apollonian packing is obtained at a ^ oo 
limit. We introduce the entropy production rate as a quantitative criteria of regularity and observe 
a transition from an irregular regime of the pattern formation to a regular one. We develop a scaling 
theory that relates kinetic and structural properties of the system. 
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A variety of physical, chemical, and biological prob- 
lems can be modeled by random sequential addition 
(RS A) processes . Examples include adhesion of pro- 
teins and colloidal particles onto surfaces, polymer chain 
reactions, car parking, etc. The structure of disordered 
media, e.g. non-crystalline and granular materials, is also 
studied within the RS A-type model approaches (cf . |^,^ ) . 
In the simplest form, RSA processes can be formulated 
as sequential addition of objects that cannot overlap and, 
once inserted, cannot move or leave the structure. 

Geometric and kinetic characteristics of RSA processes 
with identical objects are fairly well-known In con- 

trast, adsorption of mixtures has been addressed in a 
very few studies If a mixture contains a small 

number of different sizes, geometric and kinetic charac- 
teristics are primarily determined by the smallest size. In 
some applications, e.g. for adsorption of latex spheres, 
or in modeling the structures of the coal grinds, ceramic 
or metallic powders, the size distribution may spread for 
several decades . Therefore, before the smallest size 
will finally win, an interesting intermediate asymptotics 
arises. To address this intermediate regime, we consider 
a continuous distribution with sizes ranging from zero up 
to some maximal size which we set equal to unity. For 
definiteness, we choose a power-law form 



P{R) 







for i? < 1, 
for R>1, 



(1) 



for the rate of adsorption of particles with radii in the 
interval (i?, R + dR) (a > due to the normalization re- 
quirement). In the present Letter we show that RSA of 
particles with the size distribution (|l]) gives rise to frac- 
tal structures with a fractal dimension that depends on 
the exponent a. We also show that in the large a limit, 
highly regular structures arise which tend to the famous 
Apollonian packing as a — > oo. 

Let ^{t) is the fraction of uncovered area at time t 
and \E'(i?, t) is the probability that a disc of radius R can 
be placed onto the surface, with a center at some arbi- 



trary point. The definition of <f>(i) and ^'(i?,t) implies 
*(0,t) = $(t) and 



d<P 
'dt 



dRP{R)'^{R,t)ndR'^. 



(2) 



Eq. (||) is written for the general case of adsorption onto 
a d-dimensional substrate, and Qd denotes the volume of 
the d-dimensional unit ball. 

For the size distribution (|l|) one can assume a scaling 
behavior for the function ^'(i?,i), i.e. 



(3) 



Here S{t) ^ t~" is a typical gap between neighboring 
adsorbed particles, and the scaling description should 
be relevant in the scaling regime, t ^ 1 and i? <C 1 
with R/S{t) finite. Without loss of generality we can set 
F(0) = 1 and thus we get $(i) = S'^{t). Therefore, the 
scaling assumption for ^'(i?,t) implies ^{t) ^ with 
z = 6v. Expressing Eq. (^ in the scaling form yields 
V = [a + d)~^ and z ~ aVld Jq dx x°''^'''~^ F{x). Hence 
for d = 1 we have v = {1 ~\- a)^^ in agreement with the 
exact result of Ref. where it was justified that for 
the power-law distribution ([T]) , the scaling behavior does 
appear in ID. 

Now we relate the geometric (fractal) properties of the 
arising patterns with the kinetics. Let n{R) be a number 
of adsorbed particles per unit area, with radii ranging 
from Rto R + dR. Clearly, 



n{R) = / dtP{R)^{R,t) - R 



a-l+(z-l)/i^ 



(4) 



We can determine the fractal dimension Dfoi the pore 
space by introducing a cutoff size e, and calculating the 
number of particles per unit volume, N(e), with radii 
greater than e. When this number behaves as a power 
law at small sizes limit, i.e., 
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N{e) = J dRn{R) - e''^' (5) 

one concludes that the fractal dimension is Df (see e.g. 
|p2[). Notice that the porosity, i.e. the fraction of uncov- 
ered area, behaves as e'^~^f . Combining Eqs. and (||), 
we arrive to the following result for the fractal dimension: 

Df = d-z{d + a) (6) 

Eq. (|^) shows that the fractal dimension of arising pat- 
terns is intimately related to the exponent z, which de- 
scribes the pattern formation kinetics. Similar relations 
between kinetics exponents and spatial properties of the 
patterns were observed for other pattern formation mod- 
els The exact result for z in ID @ completely 
solves the problem for the one-dimcnsional systems. 

For d > 2 we could neither find the exact value of the 
exponent z, nor strictly prove the existence of scaling. 
To check the validity of the scaling theory, we performed 
numerical study of the 2D adsorption process. For every 
value of a we generated 10^ discs, and then computed the 
fractal dimension from the relation N{e) ~ e'^^f p2|-p^ . 
Typical patterns for small and large values of a are shown 
on Fig. 1 and Fig. 2, respectively. Fig. 3 gives the frac- 
tal dimension as a function of the exponent a. We also 
present the value of fractal dimension, calculated from 
the scaling relation of Eq. (|^), where numerical results 
for the kinetic exponent z were used. Fig. 3 indicates 
that the scaling theory works fairly well in the whole 
range of a. 

The lack of evident spatial correlations at small a (see 
Fig. 1) suggests that a mean-field theory (MFT) can pro- 
vide a reasonable approximation in this region. We de- 
velop a MFT similar to the one of Refs. and de- 
rived the following equation for the scaling function F{x): 

F(i?r) = (7) 

exp {^-J^dT dpP{p)a^ {{R + pY - P'') Fipr'^)^ . 

The ansatz F{x) = exp(— Aio; — ... — Adx'^) solves 
Eq. (|^) in arbitrary dimension. In particular, in 2D the 
coefRcients Ai and A2 are determined from equations 

POO 

Ai=2Tr{a + 2)a / x"'e-^''=~'^^'=^ dx 
Jo 

A2 = ^{a + 2)a / x''-^e~'^'^-'^'^' dx (8) 
2 Jo 

Solving Eqs. (^ and then inserting F{x) = ^-^i^-^^x^ 
into expression z — an dxx"~^^F{x), derived previ- 
ously, one finds z and Df. In particular, Df^2 — na — 
(tt'^/^ -I- 7r/2)a^ -I- ... in the small a limit. Mean-field re- 
sults are given in Fig. 3. Fig. 3 indicates that the MFT 
works quite well for small a, and it can be shown that 
MFT becomes exact at a ^ limit. It is failed, however, 
for large a, where the spatial correlations seem to be very 



important and the arising patterns strongly resemble the 
regular structures (see Fig. 2). 

To quantify the increasing regularity of the structures, 
we introduce an entropy, Sn, characterizing the degree 
of order of A^-particles patterns. If Ck{N) denotes such 
a pattern and p{Ck) denotes the probability of that pat- 
tern, one can define the Shannon entropy 

Sn ^ -Y.P(Ck)\og2p{Ck) (9) 

As it follows from the Eq. (||), Sn = for a regular 
pattern, since only one definite configuration with the 
probability p — I contributes to the entropy. Sn rapidly 
increases with increasing number of possible configura- 
tions, i.e. with decreasing order. The closely related 
value, dSN/dN ~ Sn+i — Sn gives the entropy produc- 
tion rate and characterizes the regularity of the pattern 
formation process. 

To find the entropy production rate, we first compute 
the conditional entropy SN[Ck{N — 1)] of the N discs 
for the given pattern Ck{N — 1) of — 1 discs. The 
conditional entropy is determined and calculated as fol- 
lows: The A^th disc added to the pattern Ck{N — 1) 
can be treated as a point in the configuration (x, y, R)- 
space, where x and y are coordinates of the disc center, 
R the disc radius. We divided the continuous configura- 
tion space into sufficiently small discrete cells and enu- 
merated these cells (for computations, we used ~ 10^ 
cells). Then, we introduced the probabilities pi that 
the next disc comes to the zth cell for the given con- 
figuration Ck{N — 1), and calculated the conditional en- 
tropy as SN[Ck{N - 1)] = -Z^iPilogjPi. The proba- 
bilities Pi can be determined numerically for any given 
configuration Ck{N — 1) of A^ — 1 discs. Namely, pi = 
A^i^R";^^ / J2j where Ai=0 if the disc correspond- 

ing to the zth cell overlaps with some disk in the pattern 
Cfc(A^ - 1), otherwise A^ = 1. 

The definitions of the conditional and full entropy then 
allows us to find the entropy production rate, Sn+i — Sn, 
by performing the averaging of the conditional entropy 
SN[Ck{N ~ 1)] over all possible configurations Cfe(A^ — 1) 
of A^ — 1 disks. In practice, the averaging was performed 
over a number of 10^) Monte-Carlo runs. 

To compare the entropy production rate for different 
values of a, we plot dS n /dN as a function of the free area 
$ (see Fig. 4). Fig. 4 shows the striking behavior of the 
entropy production rate at very large a: At the begin- 
ning of the process of pattern formation (i.e. at $ « 1), 
it decreases slowly in the same manner as for small a, 
but at <& « 0.55 a sharp decay to a (small) plateau value 
is observed. One can interpret such type of behavior as a 
transition from a regime of "low regularity" to a regime of 
"high regularity" in the pattern formation process. The 
threshold value, $ ~ 0.55, is close to the jamming den- 
sity, $00 = 0.542 . . ., of the ordinary RSA of identical 
discs on a plane. The transition can be understood as 
follows: At the first "non-regular" stage of the process. 
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random patterns mainly composed of discs with R = 1 
arise. After the jamming limit is achieved, the rule "in- 
sert the disc of maximal possible radius" starts to work 
[as it follows from Eq. (Q) in the limit a — > oo] , and much 
more regular patterns are produced. This transition cor- 
responds to the sharp decrease of the entropy production 
rate observed in Fig. 4. 

Now we argue that the fractal dimension of the pat- 
terns arising at the limit a ^ oo tends to that of the fa- 
mous Apollonian packing, the oldest known fractal con- 
struction (see To realize this we notice that in 
accordance with the rule "insert the disc of the maxi- 
mal size", that works at $ < $00, the every new added 
disc touches (at least) three other ones. Therefore, after 
some time, all the free area will be scattered into a fi- 
nite number (say L) of curvilinear triangles, confined by 
discs of radii (a;, bi, ci), / = 1, . . . , L. Contribution to the 
free area from the other curvilinear polygons (confined 
by 4, 5, etc., touching discs) is negligible if time is large 
enough. Filling of these uncovered curvilinear triangles 
is performed in the very same way as in the Apollonian 
packing, namely, the every next disc placed touches from 
the inside the three outer touching discs. The fractal di- 
mension of the pore space, obtained by the Apollonian 
filling does not depend on the starting set (a;, bi,ci) and 
equals to Da = 1.305 ... |17|. Thus we conclude that 
the fractal dimension of the patterns arising at a — > 00 
coincides with that of the Apollonian packing. 

In summary, we have investigated the adsorption ki- 
netics and spatial properties of the arising patterns in a 
random sequential adsorption of particles with power-law 
size distribution. We have developed a scaling approach 
and verified scaling description by comparing with ex- 
act results in one dimension and numerical results in 
two dimensions. We have found that arising patterns 
have the fractal dimension Df that sensitively depends 
on the power-law exponent of the particles size distri- 
bution and is intimately related to the kinetic exponent 
of the time-dependent coverage. We have deduced that 
when a, increases from to cxd, Df decreases from 2 
to Da = 1.305 . . ., of the Apollonian packing. We have 
observed that the regularity of patterns increases with 
increasing a and introduce the entropy production rate 
as a quantitative criteria of the regularity. 
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Fig. 1 A typical pattern for a = 0.1. Only a small part 
of the total number of discs is shown. 

Fig. 2 The same as in Fig. 1, a = 50. 

Fig. 3 The fractal dimension Df versus a. In the inset: 
kinetic exponent z versus a. The mean- field results are 
obtained from numerical soluion of Eqs. (||). 

Fig. 4 The entropy production rate versus the uncov- 
ered area $ for different values of a. 
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